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Statistical mechanics of bent twisted ribbons 

Supurna Sinha^ and Joseph Samuel"'^ 
^ Raman Research Institute, Bangalore, India 560 080 

We present an analytical study of bent twisted ribbons. We first describe the elastic response of 
a ribbon within a purely mechanical framework. We then study the role of thermal fluctuations in 
modifying its elastic response. We predict the moment angle relation of bent and twisted ribbons. 
Such a study is expected to shed light on the role of twist in the "J factor" for DNA looping and 
on bending elasticity of twisted graphene ribbons. 

PACS numbers: 82.35.Pq,82.37.Rs,87.10.Pq,05.20.-y 



DNA is a twist storing polymer. Many biological pro- 
cessess in the cell exert torsional stresses on the DNA 
molecule. One would expect that such torsional stresses 
affect the bending elasticity of the molecule and its abil- 
ity to form loops. "Cyclization" of DNA has been the 
focus of several studies, theoretical and experimental, 
in vitro and in vivo\^:\^. A proper understanding of 
DNA elasticity over a range of length scales is valuable 
to understanding biological processes. The persistence 
length Lp of DNA is about 50nm. At length scales 
short compared to Lp the molecule is well described by 
a purely energetic treatment similar to that used by civil 
engineers to describe the twisting and bending of beams 
and cables [5|-ll|. At length scales of the order of the 
persistence length, the effect of thermal fluctuations be- 
comes appreciable. The subject of this study is the role 
of thermal fluctuations in shaping the elastic proper- 
ties of torsionally stressed macromolecular beams. Our 
study is general and applies equally to other systems 
like actin fllaments or carbon nanotubes under mechan- 
ical stresses [8J. We use an approximation scheme [9l4ll 
developed earlier to account for thermal fluctuations 
using the Van Vleck correction. Such an approxima- 
tion works well in the stiff limit and as described ear- 
lier, works [lOj suprisingly well even over length scales 
as large as flve times the persistence length. We work 
within the wormlikc chain model 12, ll3| , which views 
the molecule as a ribbon, carrying an energy cost for 
bending and twisting. The wormlike chain model has 
been known to describe double stranded DNA [ij] as 
well as actin fllaments jl5|. Fig. 1 shows typical config- 
urations of a twisted bent ribbon. 

For clarity we consider experiments in which one end 
of the molecule is tethered to a glass slide, its tangent 
vector at the same end is constrained to lie along the 
ti direction and the molecule is torsionally constrained 
at both ends. We can attach a magnetic bead to the 



FIG. 1. Figure shows two configurations of a twisted and 
bent molecule with fixed end tangent vectors ti and tf. 
On the left (Fig. la) is the minimum energy configuration 
and on the right (Fig. lb) the same configuration is shown 
slightly perturbed by thermal fiuctuations. 



untethered end and apply bending moments and twist 
up the molecule by varying the direction of an applied 
magnetic field. Note that we do not constrain the final 
position of the molecule x{L) but only its final tangent 
vector tf . A uniform magnetic field will result in a pure 
bending moment without applying any stretching force. 
Plotting the bending moment vs the bending angle gives 
another experimental probe of the elastic properties. In 
this paper we derive the predictions of the wormlike 
chain for such experimental situations. 

We present an approximate analytical solution of 
this statistical mechanical problem. Working with the 
wormlike chain model, we first write down a closed 
form analytical solution in parametric form, express- 
ing the mechanical energy functional of a twisted bent 
ribbon as a function of the bending angle. We use this 
to plot the theoretically predicted moment-angle rela- 
tion in the absence of thermal fluctuations. We then 



compute the fluctuation determinant which describes 
the leading thermal correction to the energy functional. 
The fluctuation determinant is used to discuss the sta- 
bility of the classical solutions. Finally, we describe how 
thermal fluctuations modify the moment angle relations 
and conclude with a discussion. 

In the wormlike chain, we model the polymer by a rib- 
bon, a framed space curve {a;(s), e*(s)}. x{s) describes 
the curve, i{s) = ^, its tangent vector and e*(.s) the 
framing (the e* s are an orthonormal frame with e"^ = i) • 
s is the arc length parameter along the curve ranging 
from to L, the contour length of the curve. x(0) = 
since one end is fixed at the origin. The tangent vectors 
at both ends i{0),i{L) are fixed to ti and if respectively. 

By simple transformations described elsewhere ([ll|, 
16|,|l7[) we can reduce this problem by eliminating the 
twist degree of freedom. The reduced problem deals 
with a space curve x{s) with tangent vector t{s) subject 
to a writhe constraint. The energy functional is 



'^''>-i£'t.'>-'"'^' 



(1) 



where A is an elastic constant with dimensions of energy 
times length. The quantity W is the writhe which is 
constrained to take a fixed value. The constraint is 
enforced by the Lagrange multiplier t, which has the 
physical interpretation of the torque imposed on the 
ribbon. (The writhe of an open space curve whose initial 
and final tangent vectors are fixed is defined [17| by 
extending the curve beyond its ends to infinity adding 
straight line segments with constant tangent vectors ti 
and if.) We will work in the constant torque ensemble, 
where r is held fixed. The quantity Lp = A/kT is the 
persistence length. The mathematical problem we face 
is to compute the partition function 



Q{ii,if) = ^exp- 



£{C) 



kuT 



(2) 



In Eq.(I2]), the sum is over all allowed configurations 
of the polymer, those which satisfy the boundary con- 
ditions for the tangent vector at the two ends: ii = 

iio)jf = iiL). 

In the stiff limit, one can neglect thermal fiuctuations 
and simply minimize the energy functional appearing 
in ([2|). The problem is formally similar to that of a 
symmetric top, (as was known to Kirchoff|l8|) and this 
greatly aids the solution. Borrowing from classical me- 
chanics, we can use variational techniques and identify 
"constants of the motion" The Euler-Lagrange (E-L) 
equations that one arrives at from the variational prob- 



lem are|17|: 



i = -T(i X i) + 7i, 



(3) 



where the term jt arises since 6t ■ t = 0. The total 
"energy" 



n = ^{e^ + sm"^ 0(j>'^) 



(4) 



is a "constant of the motion" as is the "angular momen- 
tum" 



J ^ {ix i) -ri. 



(5) 



For positive (negative) r, we choose the —z {z) axis 
along the direction of the angular momentum J . Below 
we restrict to r positive. Similar considerations apply 
to negative r with appropriate changes. We have 



J=Jz 



T COS 6, 



(6) 



where we have introduced the usual polar coor- 
dinates on the space of tangent vectors i ~ 
(sin cos (/), sin 6* sin 0, cos 6*) . Evaluating —i.J = t — 
J cos 9 we find that the angle is a constant of 
the motion and that (j) ~ —J. fJ- = cos 6 is pos- 
itive (for T > 0, a condition we assume). The 
solution of the Euler-Lagrange equation is tci{s) = 
(sin 6 cos Js, — sin 6 sin Js, cos 6) . The tangent vector is 
restricted to a fixed latitude on the unit sphere (corre- 
sponding to "precession without nutation" in the top 
analogy). In real space the molecule describes a helix. 
By choice of the x axis we set the initial and fi- 
nal tangent vectors to: ti = (sin6',0, cos6') and tf = 
{sinO cos (f>f, sin 9 sin (j)f, cos 6). Considering the dot 
product u = ti.tf = cos a and setting fi = cos 9 gives us 
u and therefore a as a function of /i. u(/i) is plotted in 
Fig. 2. 



tL 

u{fi) = (1 - fi^) cos( ) + /i^ 

and an expression for a as a function of /i. 

tL 
q;(/l() = arccos [(1 — /x^) cos( — ) + fi^ 

1^ 



(7) 



(8) 



For fixed a, u = cos a defines a horizontal straight 
line (shown in black in Fig. 2), which intersects the 
curve (Eq. ^) at multiple valines of /i. The largest of 
these corresponds to the lowest value of J. The smaller 
values of ^ correspond to higher J values and higher 
rates of traversal, including traversing the circle of lat- 
itude more than once. As we will see below, only the 
largest fj. value gives a stable configuration. 
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FIG. 2. A plot of u{^) versus /i. The black horizontal line is 
a line of fixed a. The red parabolic line represents 2^^ — 1 
which is the envelope of the curve u{a) near its minima. 
Note that the minima are very close to the points of tangency 
between the envelope and u{a). 



Note that tci remains in the northern hemisphere 
since jjl > Q. Such curves (and their small perturba- 
tions) arc "good curves" [17|, in the sense that they 
can be continuously deformed to the straight line along 
z without self intersection or pointing towards —z. 
We can therefore use Fuller's formula for the writhe 
27r>V = /„ d(t>{\ — cos 6') and compute the energy func- 
tional in terms of ^ by doing the integrations in ([1]) for 
the classical solution: 



^c/(/i) = ^[(^ + l)'-4] 



(9) 



Equations (JH) and ([9|) give the solution £ci{(y) in para- 
metric form. Differentiating £ci{oi) with respect to a 
gives us the moment-angle relation: 



M^i{a) 



da 



(10) 



where the subscript d on the left hand side of (fTO| re- 
minds us that this expression ignores thermal flucua- 
tions. 

So far, we have dealt with the classical solutions of 
the Euler-Lagrange equations, ignoring questions of sta- 
bility and thermal fluctuations. Both of these ques- 
tions will now be addressed. As explained in |ll| in 
great detail, the fluctuation determinant A can be cal- 
culated from the Energy functional £ci(ti,tf) using the 
Van Vleck determinant 



A~^ = det 



d^Sc. 



dUdt 



f 



(11) 



This computation of this two by two determinant yields 
an elegant formula for the inverse of the fluctuation de- 



A-i = 



_2 d_ 

sin a da 



[Ml(a)] 



(12) 



where Mci{a) the classical Moment-angle relation. Eq. 
(|12p is the main result of this paper can be used to un- 
derstand both the stability and the fluctuations around 
the classical solutions. Eq. ([T2|) can be used to graphi- 
cally display the free energy and moment angle relation 
of the ribbon. 

The fluctuation determinant is the product of the 
eigenvalues of the fluctuation operator. A stable con- 
figuration is a local minimum of the energy and so the 
fluctuation operator has only positive eigenvalues. In- 
stability sets in when at least one of the eigenvalues 
crosses zero. At such points the fluctuation determi- 
nant A vanishes. Thus, A Eq. ([T2|) can be used to 
study the stability of classical solutions. Note first that 
H = 1 describe a straight line configuration. These con- 
figurations are known to be stable for t < tt/L [6|,|17[, a 
condition we will assume hereafter. Suppose < /i < 1, 
so that a ^ 0, TT. At points fiQ where du/dfi{fio) ^ 0, 
Taylor expansion reveals that u{p,),£tfi) and a{fi) vary 
linearly in (/i — /xq). A small calculation shows that 
A^^ is finite and so A does not vanish. It follows that 
configurations remain stable as long as du/dii{pio) ^ 0. 
Instability sets in when 



du(ii) 
dfi 



= 



(13) 



From Fig. 2 we see that the green region is stable but 
the remaining configurations arc not. Eq. (|13p is a 
transcendental equation which can be solved numeri- 
cally. A glance at Fig. 2 reveals that the roots of this 
equation are close to (but not equal to) the points where 
u(/i) touches the envelope (shown in red online) 2^^ — 1, 
where a = 29. For a ribbon bent by angle a and sub- 
ject to a torque r, the critical line for buckling is close 
to r = tt/L cos a/2, the exact value being given by the 
roots of Eq. ([T3|) . For torques and bending angles above 
the critical line the molecule buckles into a loop and cy- 
clizcs. It is evident from this formula that higher torque 
implies lower a for buckling. Thus torque aids cycliza- 
tion as could be naively expected. 

Another factor that aids cyclization is thermal fluctu- 
ations. The thermal fluctuations are easily incorporated 
in using Eq. (|12p in the expression for the free energy 

(14) 



^(a(M)) = £,i{fi) + l/2kT\ogA{fi) 



This free energy can be differentiated to give the mo- 
ment angle relations when thermal fluctuations are in- 
corporated. The results are best seen in graphical form. 



Fig. 3 shows the Hclmholtz energy as a function of a 
and Fig. 4 shows the Moment- Angle relations with and 
without thermal fluctuations. 
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ing angle and then for large bending angles it diverges. 
Correspondingly, the moment is initially Hookean and 
then diverges for large values of the bending angle. 

In the presence of thermal fluctuations there is a soft- 
ening of the "bending rigidity" which is reflected in the 
curves depicting the free energy and the moment an- 
gle relations. This has implications for cyclization and 
DNA looping "J factors". Our analysis suggests that 
in a cellular environment, torsional stresses and ther- 
mal fluctuations enhance cyclization probability of bent 
polymers. 

To summarize, this paper offers an analytical treat- 
ment of the mechanics and thermal fluctuations of 
twisted and bent polymers. We give a parametric solu- 
tion of the problem, which can be used to generate plots 
of experimentally accessible quantities. We expect that 
this solvable model would be of interest to physicists 
and biologists interested in the cylclization of DNA. 



FIG. 3. The Helmholtz free energy J-{a) versus a for 
parameters t — .3, L — 3, kT = .5, with (thick red curve) 
and witliout thermal fluctuations (thin blue curve). 
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FIG. 4. Bending moment versus bending angle for the same 
parameter values as in Fig. 3 caption (r — .3, L = 3, kT = 
.5). Thin blue line shows the mechanical elastic response and 
the thick red line includes the efi'ect of thermal fluctuations. 

To conclude, in this paper we have studied the elas- 
tic response of a bent twisted ribbon and have analyzed 
how thermal fluctuations modify its response. Note that 
the r — )■ limit is a singular one and has to be taken 
with care. From t = J^ it follows that r — ?> also 
implies /i — ;• 0. Both r and /i need to be taken to in 
such a way as to keep J finite. We find that classically 
the energy is quadratic for small values of the bend- 
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